Spike Onset Dynamics and Response Speed in Neuronal Populations 
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Recent studies of cortical neurons driven by fluctuating currents revealed cutoff frequencies for 
action potential encoding of several hundred Hz. Theoretical studies of biophysical neuron models 
have predicted a much lower cutoff frequency of the order of average firing rate or the inverse 
membrane time constant. The biophysical origin of the observed high cutoff frequencies is thus not 
well understood. Here we introduce a neuron model with dynamical action potential generation, in 
which the linear response can be analytically calculated for uncorrelated synaptic noise. We find 
that the cutoff frequencies increase to very large values when the time scale of action potential 
initiation becomes short. 
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In the cerebral cortex of the brain information is en- 
coded in the action potential (AP) firing rates of a large 
ensemble of nerve cells. Recent experiments have ob- 
served a surprisingly high cutoff frequency for the action 
potential encoding of cortical neurons driven by fluctuat- 
ing input currents [H-0| ■ I n a seminal paper Kondgen et 
al. showed that the transmission function of layer 5 pyra- 
midal neurons for a noisy sinusoidal signal does not de- 
cay until about 200 Hz [I| . Later experiments confirmed 
such high cutoff frequencies for signals coded by both the 
mean current and noise strength Q and in other types 
of cortical neurons [3, 4]. For an early observation of fast 
response see [B| . Previous theoretical studies of biophysi- 
cal neuron models, however, predicted cutoff frequencies 
of the order of the average firing rate or the inverse mem- 
brane time constant (below 20 Hz), much lower than the 
experimentally observed values [6|-l8|- Thus, the origin of 
the high cutoff frequencies found in cortical neurons is 
currently not well understood. Numerical investigation 
of neuron models with dynamical AP generation, like the 
exponential integrate-and-fire (EIF) model or the gener- 
alized theta neurons, suggested that details of AP gen- 
eration can influence the dynamical response of neuronal 
populations What is missing, however, is a trans- 

parent understanding of how and when the population 
cutoff frequency can dissociate from the basic single neu- 
ron timescale set by the mean firing rate and the time 
constant of membrane potential relaxation. 

In this work we present an analytically solvable model 
which explicitly describes the dynamical AP initiation 
process. A neuron initiates an AP if the membrane po- 
tential passes an unstable fixed point, the voltage thresh- 
old. In the leaky integrate-and-fire (LIF) model, for 
which the linear response is known analytically, the un- 
stable fixed point coincides with the absorbing bound- 
ary and a spike is triggered immediately when the mem- 
brane potential reaches this threshold [l(| El • As a con- 
sequence, boundary induced artifacts dominate the_ re- 
sponse for high signal frequencies in the LIF model 



One important advantage of our new model is that such 
boundary induced artifacts can be separated out mathe- 
matically, isolating the physically meaningful part of the 
response function. We first present the linear response 
for both encoding paradigms with white noise. We find 
that for a wide range of parameter settings the cutoff 
frequency is directly proportional to the AP onset rapid- 
ness for a noise coded signal. It therefore dissociates from 
the membrane time constant and can become arbitrarily 
large. For the mean current coded signal, however, the 
cutoff frequency is confined by the membrane time con- 
stant in the white noise case. We show by numerical 
simulation that this confinement can be broken when a 
finite correlation time in the synaptic noise is taken into 
account and high cutoff frequencies can be obtained for 
a large AP onset rapidness. Interestingly, experiments 
showed that the AP onset rapidness of cortical neurons 
is very large both in vitro and in vivo [l2|j EH; which 
may thus explain the occurrence of high cutoff frequen- 
cies. Our results provide a relationship between the spike 
onset dynamics and the population cutoff frequency that 
can be directly tested in physiological experiments. 

The simplest voltage dynamics that exhibits both a 
stable fixed point (the resting potential) and an unsta- 
ble fixed point (the voltage threshold) has a piecewise 
linear membrane current, composed of a leak current for 
low potential and a linear spike generating current for 
high potential [Fig. QJa)]. The model is defined by the 
following Langevin equation 

r m v = -v + Q(v - v a ) (r + l){v - Vq) + [1 + GT)(t) , (1) 

where v is the membrane potential relative to the rest- 
ing potential, r m is the membrane time constant, ®(v) 
is the Heaviside step function, r is the AP onset rapid- 
ness which sets an effective time constant r m /r for the 
AP initiation process. The larger is r, the faster is the 
spike onset [Fig. [IJb)]. In biophysical models, the on- 
set rapidness will be set largely by intrinsic properties of 
the voltage dependent sodium channels, e.g., the gating 
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charge and the slope of the activation curve 12| . \i is the 
mean input current and a is the amplitude of synaptic 
noise. 77(f) is a Gaussian white noise satisfying (77(f)) = 
and {r)(t)r](t )) = T m S(t — f ). The crossing point i> of 
the two pieces sets the rheobase current, which we use as 
the unit of voltage, Vq = 1. The threshold potential is 
vt = (1 + l/r)VQ. When the membrane potential reaches 
Vb, the truncation point of the AP upstroke, it is reset 
to a voltage v r and stays there for an absolute refractory 
period t t . For convenience we take r m as the unit of time 
in analytical calculation. 

The Fokker-Plank equation corresponding to Eq. (JTJ 
has the following form 

d t Pi +dj-v + n- i^a^Pi = , 

d t P 2 + d v (r(v -v t )+ii- \cx 2 d^j P 2 = 0, (2) 

where Pi(v, f) and P%{v, f) are the probability densities 
of membrane potential v for — 00 < v < 1 and 1 < v < Vb 
respectively. The stationary firing rate can be found eas- 
ily when the boundary conditions are specified. From the 
reset assumption, we impose an absorbing boundary at 
Vb and, as a result the probability density is zero there, 
Pzivbjt) — 0. The firing rate is given by the probabil- 
ity current at w&, v(t) = — ^a 2 d v P 2 (vb,t). At the reset 
voltage v r , the probability density is continuous, while 
its first derivative has a discontinuity from the reset con- 
dition: d v Pi(v+,t) — d v Pi(v~,t) = d v P 2 (vb,t — r r ). In 
addition, the density and its first derivative should be 
continuous at v = 1. When /j and a are constants, the 
system is homogeneous and the stationary solution of 
Eq. ((2]), denoted as Poi(v) and Pq2(v) respectively, can 
be found. The stationary firing rate vo is then obtained 
from the normalization condition of the density 

The stationary firing rate vq of the model Eq. 
([I]) reduces to that for the LIF model for r — > 00. Figure 
[T](c) and D show the dependence of u on the mean input 
(i and the amplitude of noise a, respectively. The firing 
rate Uq increases monotonically with r, //, and a and is 
relatively insensitive to r when r > 10. For the dynam- 
ical response, however, the r dependence is much more 
pronounced. The instantaneous firing rate of an ensem- 
ble of model neurons responds much faster for larger r to 
a step change in the noise level 

Linear response. — When the input current to a neu- 
ron is weakly modulated, linear response theory can be 
applied to study the dynamical response properties of an 
ensemble of neurons. To this end, we consider a sinu- 
soidal signal ecos(wf), where e is small. When the signal 
is encoded in the mean current, fj, — > /j + ecos(wf) or in 
the noise amplitude, a —> a + ecos(o;f) , the instantaneous 
firing rate can be written as v(t) = i'o-l-e|i / i c (w)|cos(a;f — 
<fi c (w)) or i/(t) = uq + ecr|^i I i(a;)|cos(a;f — <p n (u:)). Here 
v\(lo) are the complex response functions. The absolute 
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FIG. 1: (color online) (a) Illustration of the model, (b) V(t) 
trajectories for identical noise and three different values of 
r. (c) and (d) show the dependence of stationary firing rate 
on mean input current and noise strength in the noise driven 
regime. (r m = 10 ms, v r = 0, v b — 10, r r = 0) 



value \ v\ (cj) | are the transmission functions and the phase 
angles 4>i(lo) — arg(z/i(a;)) give the phase lags, which 
completely characterize the linear response. Note that 
we refer to both signal channels when the subscripts c 
and n are omitted. 

It is known that the absorbing boundary condition 
at Vb can induce severe artifacts in the dynamical re- 
sponse. The potential Vb marks a "point of no re- 
turn," which is not present in a biophysical dynamical 
model of AP initiation. As a consequence, the trans- 
mission function for a noise coded signal in the LIF 
model, for instance, does not decay at high signal fre- 
quencies [H [ll|. Ideally one would thus wish to sepa- 
rate the response function into a physiologically mean- 
ingful part v (cj) and a part containing all artifacts 
such that vf hy (lu) = v\(lo) — z/° bs (w). v^ hy (oj) must have 
the following properties: i) vf v (oj) approaches the static 
susceptibility when uj — > 0, specifically, i/^ y (u>) — > 



and <7(w) -> if^; ii) ^ ny (uj) -> when u -> 00; 
iii) no essential dependence on the truncation point Vb- 
The artifactual part from the absorbing boundary should 
have the following properties: i) negligible contribution 
for signal frequency in the physiologically relevant range 
/ < 1kHz, e.g. \vf s {io)\ < \v{ hy {oj)l where / = w/2tt; 
ii) strong dependence on the truncation point Vb- As we 
will show next such an isolation of the physiologically 
meaningful response is possible in the model Eq. <jTJ) . 

In the model Eq. (fTJ, the linear response can be ob- 
tained analytically by expanding the probability density 
in Eq. ([2]) to first order in e and using the Green's func- 
tion method. We find that v\ n (ui) decomposes naturally 
into two parts, v\(uj) = i>i° w (u]) + v^ tgh {uj), with 
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tributions imposed by the absorbing boundary condi- 
tion is finally confirmed from the high frequency behav- 
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quency behavior in the LIF model is determined solely 
by the absorbing boundary. As a consequence, neglect- 
ing v^J 3 (u)) in the response function eliminates any 
boundary induced instantaneous response components. 
These results establish that v^ ow (lo) capture the be- 
havior of V\ (iS) for low and intermediate frequencies 
and decay to zero in the large frequency limit. When 
M -> v ln (0) = i^, since i£* h ( W ) is 
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FIG. 2: (color online) The normalized function v\(iS) jv\{f).V) 
and phase lag for a mean coded signal and noise coded signal 
with r = 1, fi — 0, and = 5 Hz. Other parameters are the 
same as in Fig. [T] Lines, theory; Symbols, simulation. 



negligible there. v\ ow (yi) exhibits only a weak depen- 
dence on Vb through a frequency dependence phase lag 
4>o — 7 l°g V ^/ + d > characterizing the time lag to the trun- 
cation point vi, of the AP upstroke. Therefore, we have 
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(lo). The physiologically meaningful pre- 



dictions of the model can thus be revealed by examining 



u^ ow (u}) in isolation. 



L Cutoff frequency and AP onset rapidness. — Figure [3] 

where D = ^a 2 , A = (1 - v r )(2fi - 1 - v r )/4D and shows the behavior of vf-^foj) with increasing r and how 
Ai = (1 — v b )(2fi - 1 + r(v b - «t))/4D. tpi, $i, and Ti 
are parabolic cylinder functions, and Yy, Y% are a combi- 
nation of parabolic cylinder functions, whose dehnition 
together with the v± ig (uj) parts are given in the sup- 
plement 14} . Any prime represents the derivative with 
respect to v. Here and in the following, the functions 
adopt their values at v — 1 if not denoted explicitly. Fig- 
ure [5] illustrates the linear response with r = 1 as an 
example. 

Removing boundary induced artifacts. — The decompo- 
sition of Vx(ui) into two additive components has exactly 
the features required for the separation of artifacts. Us- 
ing asymptotic expansion of the parabolic cylinder func- 
tions we find that for a finite signal frequency, v^ %gh can 
be approximated by 



the cutoff frequency f c changes with r for different vq 
We see that f c increases linearly with the onset 
rapidness r for noise coded signals when the firing rate 
is not very low (> 1 Hz here); while for mean coded 
signals f c saturate for large r. The increase of f c with 
firing rate v$ results from the stochastic double resonance 
phenomenon: the transmission function will develop a 
peak for some optimal signal frequency before decaying 
when vo is relatively large 16]. 

Figure [3)3 suggests that the cutoff frequency for a 
noise coded signal follows f c oc r and dissociates from 
T m . This is directly confirmed by the large frequency 
approximation of Vy° w (uS), 
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when Vb ^> v t . So v^ lgh (uj) are strongly dependent on 
Vb and approach zero when Vb — > oo for finite signal 
frequency. When u —> 0, v^ lgh {ui) are negligible com- 
pared with Vi° w (uj) and are strongly suppressed when 
Vb is large. That i^ ln * s (w) captures all artificial con- 



where Poi(v) = y/rPQi(v). Because Poi(uo 

and P 01 (uo) - 
depends essentially only on cu/r. This implies that the 
cutoff frequency for a noise coded signal dissociates from 
T m and becomes proportional to the onset rapidness r, 
f c = At, where A depends on vq through the effect dis- 
cussed above. This demonstrates that fast onset APs can 
enhance the cutoff frequency and therefore the response 
speed significantly. Note that a linear relationship f c oc r 
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FIG. 3: (color online) (a) and (c) The normalized transmis- 
sion function Ui" w (u})/vi(0.1) for different r with fi — 0, 
Uo = 5 Hz. (b) and (d) The variation of cutoff frequency 
with r for different firing rates: vq — 1,5,10,20,30,40 Hz 
from lower to upper curves. Other parameters are the same 
as in Fig. [1] 



membrane time constant. The confinement of the mean 
response for white noise constitutes an interesting predic- 
tion of the model, which should be tested experimentally 
by using a very short correlation time (< 1 ms) of the 
synaptic noise. The origin of the large onset rapidness 
seen in cortical neurons is a matter of ongoing debate 
[ItI [l8| . Its value can be modified in real neurons by 
applying dru gs li ke TTX or knockout of sodium channel 
subtypes 



12, 



19j . Measurements of dynamical response 
for cortical neurons under such manipulations are thus 
predicted to provide important insight into the mecha- 
nism of fast population coding. 
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FIG. 4: Variation of the normalized transmission function 
|^i c (tj)/^i c (0.1)| with increasing correlation time r s (with 
unit ms) of the synaptic noise for r = 10, 100. Parameter 
used: jj, = 0, vt = 10, vo = 5Hz. 



was previously conjectured by Naundorf et al. based on 
dimensional analysis 18 1. 



For a current coded signal, however, 



exp(— juj/r) for large r. Therefore the linear response 
is confined by the membrane time constant in the white 
noise case, as seen also from Fig. [3j^. Real synaptic in- 
puts, however, have a finite correlation time and can be 
approximated better as a colored noise. As shown in 
Fig. 01 the confinement by the membrane time constant 
is broken under such conditions and cutoff frequencies of 
several hundred Hz can be reached for large r 14 1. 

Our results identify AP onset rapidness as a critical de- 
terminant of population cutoff frequency and reveal how 
this cutoff frequency can dissociate from the basic single 
neuron time constants set by the mean firing rate and 
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